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ABSTRACT 

We report here on the discovery of stellar occupations, observed with Kepler, that recur periodically 
at 15.685 hour intervals, but which vary in depth from a maximum of 1.3% to a minimum that can 
be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a V = 16 magnitude 
K dwarf with T c g tS ~ 4400 K. The out-of-occultation behavior shows no evidence for ellipsoidal light 
variations, indicating that the mass of the orbiting object is less than ^3 Mj (for an orbital period 
of 15.7 hr). Because the eclipse depths are highly variable, they cannot be due solely to transits 
of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant 
planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing 
transit. This scenario seems ruled out by the dynamical instability that would result from such a 
configuration. We also briefly consider an eclipsing binary, possibly containing an accretion disk, that 
either orbits KIC 12557548 in a hierarchical triple configuration or is nearby on the sky, but we find 
such a scenario inadequate to reproduce the observations. The much more likely explanation — but 
one which still requires more quantitative development — involves macroscopic particles escaping the 
atmosphere of a slowly disintegrating planet not much larger than Mercury in size. The particles 
could take the form of micron-sized pyroxene or aluminum oxide dust grains. The planetary surface 
is hot enough to sublimate and create a high-Z atmosphere; this atmosphere may be loaded with 
dust via cloud condensation or explosive volcanism. Atmospheric gas escapes the planet via a Parker- 
type thermal wind, dragging dust grains with it. We infer a mass loss rate from the observations 
of order 1 M^/Gyr, with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1M® planet 
(twice the mass of Mercury), the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are 
disfavored because they evaporate still more quickly, as are larger mass planets because they have 
surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation 
profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; 
we present simulations of such a tail. 

Subject headings: eclipses — occultations — planets and satellites: general — stars: planetary systems 



1. INTRODUCTION 

The Kepler mission has now discovered some ~1250 
very good exoplanet candidates via planetary transits of 
the parent star (Borucki et al. 2011). The transit depths 
range from a few percent for giant gas planets to < 10 -4 
for smaller planets that are of the size of the Earth. The 
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distribution of orbital periods for these systems is shown 
in Fig. IT} The bulk of the systems have periods between 
^3 and30 days. As observations continue to be accumu- 
lated and analyzed, the newly discovered candidates will 
tend to be longer period systems; the short-period end 
of the period distribution will change relatively slowly. 
There are relatively few systems with periods below a 
few days, and furthermore, most of the systems shown 
in Fig. [TJ with orbital periods less than a day are either 
now considered to be false positives (i.e., they are most 
likely associated with background blended binaries) or 
have not yet been followed up in the detail required for 
conclusive categorization (Borucki et al. 2011). The con- 
firmed short-period planets, most of which were not dis- 
covered by Kepler] 11 ] have orbital periods between 0.73 
and 1.09 days. These orbital periods are indicated by 
dual arrows in Fig. [TJ 

The shortest planetary orbital period that has been re- 
ported is 17.8 hours (Winn et al. 2011)0 The planet 55 
Cnc e, the innermost of 5 planets revolving around 55 
Cnc, was discovered with ground-based optical Doppler 

11 Taken from the "Exoplanet Orbit Database" at 
http://exoplanets.org, produced by Jason Wright, Geoff Marcy, 
and the California Flanet Survey consortium. 

12 Only just recently an exoplanet with an even shorter orbital 
period, KOI 961.02, with P or b = 10.9 hours, was reported (Muir- 
head et al. 2012). 
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Figure 1. Distribution of orbital periods among the exoplanet 
candidates discovered with Kepler. The shorter period planets 
specifically named are not necessarily from the Kepler survey and 
are taken from the on-line catalog "The Exoplanet Orbit Database" 
(sec text). 



measurements (Fischer et al. 2008; Dawson & Fabrycky 

2010) and detected in transit by extensive satellite ob- 
servations with the MOST satellite (Winn et al. 2011). 
Winn et al. (2011) estimate its mass and radius to be 
^9 and 2 . Planet orbits considerably shorter 
than a day are, in principle, quite possible, and their 
detection via transits using Kepler should not be prob- 
lematic in a signal-to-noise sense. There are at least two 
reasons, however, why such close-in, short-period plan- 
ets have not been heretofore emphasized. First, there are 
many short-period stellar binaries that, when their light 
is highly diluted by another star within a few arc seconds 
(e.g., the target star under observation by Kepler), are 
hard to distinguish from stars with transiting planets. 
Second, giant planets are expected to be destroyed by 
Roche lobe overflow if they are heated too strongly by 
their parent stars (Batygin, Stevenson, & Bodenheimer 

2011) . 

With Kepler data, the orbital period P or b of a planet 
about a star can be measured. If the orbit is circular, 
and the planet has a mass M p much less than the mass 
of the parent star, then the size of its Roche lobe can 
be determined without knowing the mass of the parent 
star: 
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where Mj and Rj are the mass and radius of Jupiter, 
respectively. Therefore, even if the planet does not evap- 
orate because of the strong stellar radiation field, by the 
time a Jupiter-mass planet migrates to an orbit with 
-Porb ^ 9 hours, it will start to empty its envelope via 
Roche- lobe overflow. 

In this work we discuss the signature of a planet in a 
15.685-hour orbit about a K star. The object came to 



our attention via periodic occultations of the target star 
KIC 12557548. The term "occultation" is perhaps more 
appropriate than "transit" (although we will use both 
terms interchangeably in this paper), because of the fact 
that the depths of the events are highly time dependent. 
As we discuss in this work, the occultations are periodic 
to 1 part in 10 5 or more, and therefore must be due to the 
presence of an orbiting companion. However, the varia- 
tion of the eclipse depths rules out transits by a single 
opaque body. More likely, we argue, the occultations are 
due to a stream of dust particles resulting from the slow 
disintegration of the planet. 

In |2.2f[2.5| of this work we describe detailed analyses 
of the Kepler observations of KIC 12557548 from quar- 
ters Q2 through Q6. In these sections we present light 
curves as well as power density spectra of the light curves. 
We also present an optical spectrum which basically con- 
firms the spectral properties of KIC 12557548 given in 
the Kepler Input Catalog (KIC). In ^3] we outline three 
explanations for the variable occultation depths, includ- 
ing grazing transits due to a pair of planets; an eclipsing 
binary, possibly containing an accretion disk, that orbits 
KIC 12557548 in a hierarchical triple configuration; or, 
more likely, obscuration by a stream of debris from a 
disintegrating planet. We explore in <H] the hypothesis 
of an occulting dust cloud in some depth. In particular, 
we discuss how dust can be stripped off a planet with 
substantial gravity. We go on to compute the approxi- 
mate shape of a dust stream that has left the planet and 
is shaped by radiation pressure. Our findings are sum- 
marized in Sj5j there we also suggest future observations 
that may robustly establish the nature of the variable 
occultation depths. 

2. DATA ANALYSIS 
2.1. Discovery 

The occultations of KIC 12557548 were discovered as 
part of a systematic search through the Kepler data base 
for eclipsing binaries that are not in the Prsa et al. (2011) 
catalog. We used publicly available quarter 2 ("Q2") 
data from each of the ^150,000 Kepler targets to carry 
out a general Fourier transform search for periodicities 
(Kotson et al. 2012). The criteria for declaring a partic- 
ular power density spectrum (PDS) as "interesting" are 
that (i) there is at least one peak with amplitude greater 
than 5 standard deviations above the local mean of the 
spectrum, and (ii) there is at least one harmonic or sub- 
harmonic of the largest peak that is at least 3 standard 
deviations above the local mean. 

In all, we find that ^4800 power spectra contain in- 
teresting content, as defined above. Of these, approxi- 
mately 1800 are binaries that are already in the Prsa et 
al. (2011) catalog, as well as some 300 binaries, mostly of 
the contact variety, that appear to have not made it into 
that catalog (Kotson et al. 2012). Most of the remaining 
FFTs of interest involve stars that exhibit pulsations of 
one sort or another. 

One of the interesting power spectra contains a peak at 
15.6854 hours and 15 harmonics thereof with amplitudes 
that decrease slowly with frequency. This is indicative of 
a sharp transit/eclipse/occultation which occupies only 
a small fraction of a cycle. The target, KIC 12557548, 
does not appear in either the Prsa et al. (2011) binary 
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Figure 2. Light curves for KIC 12557548 spanning the Q2-Q6 data sets, as well as for 40-day and 15-day segments of the Q2 data. The 
intrinsic Poisson fluctuations per ~30-min Kepler sample are ~0.03%. 



star catalog or the Borucki et al. (2011) planet catalog. 

2.2. Light Curves 

Light curves for KIC 12557548 are shown in Fig.[2j 
The plots shown in this figure were made by combin- 
ing the raw (SAP_FLUX) long-cadence light curves from 
quarters 2 through 6 (Q2-Q6); before this was done the 
light curves from the various different quarters were mul- 
tiplied by relative scale factors chosen to minimize any 
differences in intensity across quarter boundaries. Low 
frequency variations caused by instrumental, and possi- 
bly intrinsic source, variations were then removed by con- 
volving the flux data with a boxcar of duration 0.65356 
days (15.5854 hours) and then subtracting the convolved 



flux curve from the pre-convolution light curve. This pro- 
cessing should largely preserve the actual source behavior 
on orbital and suborbital timescalesPH 

Fig. [2] shows the highly variable nature of the occulta- 
tions which, in fact, have depths that render them easily 
seen over most of the 460 day interval, but for the first 
^22 days of Q2 are very small or even not noticeable. We 
have also checked the Ql data and, indeed, the occupa- 
tions are not apparent during the first 10 days as well as 
the final 10 days of that quarter. As seen in panel (c) 

1 3 We also used the P yKE sof tware to process the data. 
See http : //keplergo . arc .nasa. gov/ContributedSof twarePyKEP . 
shtml 
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of Fig.[2j the individual occultations vary strongly from 
one cycle to the next. The variable occultations persist, 
with a similar range in depths, throughout the following 
438 days of Q2 through Q6 data. 

When the data are folded about the 15.685-hour pe- 
riod, the results, shown in Fig.[3j further illustrate how 
the occultation depths vary. Thelargest depth is ~1.3%, 
while some of the points near the middle of the occulta- 
tion are depressed by no more than ~0.15%. On the 
other hand, it is clear that, aside from a few individual 
data points within the occultation interval, almost all 
show a minimum depression in the intensity by ~ 0.1%. 
Thus, even when the occultations are not noticeable, 
there still may be a small effect. 

When the folded data are binned (see Fig. [3}}) one can 
see that the average depth at the heart of the occultation 
is only ^0.6% or about half the maximum depth. There 
is evidence (at the ~5-<r level) for a small peak in the 
flux of amplitude ~3 x 10~ 4 at phases up to ^30° before 
the occultation ingress. There is a also a very significant 
depression or deficit in the expected flux over a similar 
range of phases just after egre ss. We will discuss these 



features briefly in Section £4.7 

The full width of the occultation, with the exception 
of these small features, is 0.1 of the orbital cycle. This 
corresponds to ~1.5 hours in duration, or just 3 Kepler 
long-cadence integration times. If this duration is inter- 
preted simply as indicating the sum of the radii of the 
occulting "bodies" it corresponds to (i?i + R2)/a ~ 0.3. 
However, if we take into account approximately the effect 
of the finite integration time, then (i?i + Rj2)/a ~ 0.24. 
Note that these estimates assume, without justification, 
equatorial as opposed to grazing occultations. 

We have fitted a constant plus a cosine with a period 
of 15.685 hr and a cosine of half that period to the folded 
light curve, excluding points inside the occultation inter- 
val (i.e., at orbital phases of ±0.085 cycles around mid 
occultation), in order to search for possible ellipsoidal 
light variations. We find only an upper limit of ^5 x 10 -5 
for the amplitude of such features. This limit will be dis- 
cussed below in the context of setting a constraint on the 
mass of any body orbiting the target K star. 

Finally, we find the epoch of mid-occultation in order 
to project ahead to future observations: 



t n = 55399.824(2) + n x 0.65356(1) MJD 



(2) 



where the numbers in parentheses are the respective un- 
certainties in the last significant figure given. 

2.3. Fourier Search for Modulations 

As a check for periodic modulations of the occulta- 
tion depths for the 15.685-hour period, we carried out 
an FFT of the Q2 plus Q3 data but with the portions 
of the light curve away from the occultation set equal 
to the mean out-of-occultation intensity. The purpose of 
this latter step is to suppress the noise level, without sac- 
rificing any significant fraction of the signal. The results 
are shown in Fig.|4j All 16 harmonics of the 15.865-hour 
period, out to the Nyquist limit, are clearly visible. In 
addition, a careful inspection of the amplitudes between 
the harmonics indicates some evidence for low-amplitude 
modulation-induced sidebands. However, at least a num- 
ber of these can be reproduced in an FFT of the window 
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Figure 3. Folded light curves of KIC 12557548 about the 15.685- 
hour occultation period for the Q2-Q6 data sets. Top panel - 
unbinned data; bottom panel - folded data averaged into 96 dis- 
crete bins (3 bins ~ 1 long-cadence Kepler integration time). Short 
illustrative vertical bar on the left is ± the standard error of the 
data points within a bin. Note the highly statistically significant 
depressed flux level following the main occultation. 



function associated with the occultations. Thus, we find 
no compelling evidence for periodic modulation of the 
occultation depths. 

2.4. Checks on the Validity of the Data 

Because of the unusual exoplanet phenomenon pre- 
sented in this work, we need to be especially careful to 
ensure that there are no spurious artifacts in the Ke- 
pler data for this object. In this regard, we performed a 
number of tests on the data. 
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Figure 4. Power spectrum of the detrended flux data for KIC 
12557548 with the out-of-occultation region set equal to the mean 
flux level. See text for details. 



First, as mentioned above, we checked that the oc- 
cultations are present in all of the quarters of released 
(i.e., public) data (Q1-Q6), and that the behavior of the 
occultations does not change abruptly across quarterly 
boundaries. 

Second, we investigated whether the photometric vari- 
ations of KIC 12557548 could be from another star or 
source on the sky. We examined the Digitized Sky Sur- 
vey images and found that there are no especially bright 
stars close to KIC 12557548 whose light or transferred 
charge might introduce spurious signals into the data 
stream (i.e., none with V < 13 within 3' and none with 
V < 10 within 8'). An examination of a Kepler Full 
Frame Image for Ql indicates that there are no obvious 
stars within ~38" (~10 pixels). The Kepler Input Cata- 
log does list two Kp=19 stars located 14" and 17" from 
KIC 12557548 (approximately 3.5 and 4.25 pixels), but 
these could not be the source of the photometric varia- 
tions for two reasons, (i) We checked difference images 
for each of Q1-Q6 and the apparent source of the photo- 
metric variations is coincident with the stellar image of 
KIC 12557548 (Batalha et al. 2010). Difference images 
compare the in-transit frames with out-of-transit frames 
in the neighborhood of each transit and these are av- 
eraged across all events in each quarter, (ii) Analysis 
of the correlation of the astrometry (brightness weighted 
centroids) with the photometric signature of the occulta- 
tions indicates that the position of KIC 12557548 shifts 
by no more than ^0.2 millipixels in either row or col- 
umn during any quarter. This indicates that there is 
no other background source that is producing the oc- 
cultations and that is offset by ^2" or more from KIC 
12557548, given the depth of the events and the bright- 
ness of KIC 12557548 (Jenkins et al. 2010). 

As a further check against video crosstalk from the ad- 
jacent CCD readout channels, we examined areas on the 
sky located at the same row/column position on these 
adjacent channels, and found no stars within 11" of the 
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Figure 5. Optical spectrum of KIC 12557548 taken on December 
16, 2011 at the OMM Observatory. The spectrum is a stack of 
four 15-minute exposures through ~1.6 air masses. The spectral 
features are indicative of a mid-KV star. 



row/column position of KIC 12557548. The fact that the 
video crosstalk coefficients vary by typically 50% or more 
from quarter to quarter as the stars are rotated about the 
focal plane, and can be either negative or positive, and 
the fact that the photometric signatures of the occul- 
tations are consistent from quarter to quarter, provide 
additional confidence that it is highly unlikely that the 
photometric variations are due to a source on an adjacent 
CCD readout channel. 

Third, there are no binaries in the Prsa et al. (2011) 
catalog or planets in the Borucki et al. (2001) catalog 
with periods that match that of KIC 12557548. 

2.5. Optical Spectrum 

KIC 12557548 is a ~16th magnitude star located 
at RA(J2000) = 19 h 23 m 51.89 s and Dec(J2000) = 
51° 30' 17.0". We were able to obtain an optical spec- 
trum of it on 2011 December 16 using the Observa- 
toire Astronomique du Mont-Megantic's 1.6-m telescope. 
Four 15-minute exposures were taken through 1.4 to 1.7 
air masses soon after astronomical twilight. The spec- 
trometer utilized the STA0520 blue CCD and a 600- 
line/mm grating, yielding an effective spectral resolution 
of 4.3 A. 

The four spectra were background subtracted, cor- 
rected for the quantum efficiency of the detector and the 
reflectance of the grating, corrected for Rayleigh scat- 
tering in the atmosphere, and co-added to produce the 
spectrum shown in Fig. [5] We did not carry out absolute 
flux calibrations since the purpose of this observation was 
simply to confirm the spectral class of the Kepler target 
and to check for any unexpected spectral features. 

Despite the relatively low signal to noise ratio of 
this spectrum, taken under non-ideal conditions, we can 
clearly see the signatures of a mid-K star. Comparisons 
with standard stars (see, e.g., Jacoby, Hunter, & Chris- 
tian 1984) indicate a spectral type of ^K4V. Further- 
more, we find that the equivalent width of the Na I line is 
3.6±0.4 Ap*] which is indicative of a K5-K7 star (Jaschek 

14 This is actually an upper limit to the EW because the ISM is 
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& Jaschek 1987). From this range of spectral types we in- 



fer T, 



eff,s 



4300 ± 250 K, a mass of M. ~ 0.70 



+0.08 



-0.04 M G> 



and a radius of R s ~ 0.65+q'q| Rq, assuming that the 
star is close to the ZAMS. These are all quite consistent 
with the values listed in the KIC. 

The KIC lists magnitudes of g=16.7, r=15.6, i=15.3, 
and z= 15.1 for star 12557548. From this we estimate 
the V magnitude is 16.2 (Windhorst et al. 1991). If we 
take the absolute visual magnitude of the K star to be 
+7.6 (with L s ~ 0.14 Lq), and correct for an estimated 
Ay — 0.22, this puts the object at a distance of about 
470 pc. 

3. INTERPRETATION: GENERAL REMARKS 

The data analysis presented for KIC 12557548 in the 
previous sections reveals a highly periodic set of occulta- 
tions, likely of the 16th magnitude K-star target of this 
particular Kepler field. These occultations are highly 
variable in depth, and are therefore very different from 
anything yet reported about the other couple of thousand 
transiting exoplanets (including candidates). 

The remarkable stability of the occultation period 
strongly suggests that the occultations are due to an or- 
bital companion to KIC 12557548. In |3.1| we constrain 
the separation and mass of such a companion. These con- 
straints are quite independent of any particular model. 
We then outline more specific scenarios to explain the 
variable nature of the occultations. These are: (1) di- 
rect occultations by one planet, modulated by a second 
planet that causes a prec ession of the orbital plane of 
the occulting planet ( |3.2[ ), and (2) dust that emanates 
from a single disintegrating planet ( p.3| . We consider 
both ideas qualitatively here, and then m Qwe elaborate 
on the more likely one: that the occultations are due to 
dust emitted directly from a planetary atmosphere. In 
§3.4| we discuss the alternative possibility that the oc- 
cultations are not of the K star, but are rather due to a 
low-mass eclipsing stellar binary, possibly containing an 
accretion disk. 

3.1. Constraints on Separation and Mass 
of Orbital Companion 

The regions of the light curve folded at the 15.685- 
hour period that are away from the occultations show 
no evidence — with a limiting amplitude of ~5 parts in 
10 5 — for ellipsoidal light variations (ELVs) with period- 
icities of 15.685 hr or half that. Since the amplitudes of 
ELVs are of order ELV ~ (R s /a) 3 (M p /M s ) (where the 
subscript s refers to the parent star) we can set a limit 
on the mass of the companion M p (assumed < M s ) if we 
have an estimate of the companion's orbital semimajor 
axis a. Using the parameters for KIC 12 5575 48 taken 
from the KIC, as well as our spectrum (see { 2.5 1, we find 
T eSiS ~ 4400 K, log g ~ 4.63, and R s ~ O.6577/ . These 
properties all point strongly toward a mid-K star of mass 
M s ~ 0.7 Mq. The semimajor axis of a small body in an 
assumed P Q rb = 15.7-hour orbit about such a star is 
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0.7 M, 







1/3 



Rr. 



0.013 AU. 



(3) 



We then have a rather good estimate of Rs/a of 0.23. 

expected to make a small contribution. 



When this value and our upper limit on the ELV ampli- 
tude are used in the above expression for the ELV ampli- 
tude, we obtain a constraint on the mass of the orbiting 
companion of 

M p < 3 Mj . 

Thus, it is unlikely that the parent star has another star, 
or even a brown dwarf, orbiting it at a distance of only 
2.8 -R©rj This leaves only planetary-mass companions 
as the direct or indirect cause of the occultations. 

3.2. The Dual Planet Hypothesis 



From our analysis in S3.1 a planet similar to Jupiter 
in mass and radius could be compatible with the lack of 
ELVs. Such a giant planet would also generate maximum 
eclipse depths of ~1%, like those observed. The challenge 
lies in having such a Jovian-mass planet produce highly 
variable transit depths. 

We could try to make such a model work by supposing 
that the transit is grazing, and that precession of the or- 
bital plane alters the transit geometry so that the eclipse 
depth varies. However, orbital precession induced either 
by tidal distortions or by another planet in an indepen- 
dent orbit generally occurs over timescales much longer 
than the orbital period. By contrast, the eclipse depths of 
KIC 12557548 vary markedly from orbit to orbit. More- 
over, invoking a second planet in an attempt to force 
rapid orbit-to-orbit variations in the planet's trajectory 
across the star would result in transit timing variations 
(TTVs; Agol et al. 2005; Holman & Murray 2005). No 
TTVs are observed, although the long-cadence K epler in- 
tegration time hinders us from detecting TTVs substan- 
tially shorter than ~30 minutes. In addition, the peculiar 
ingress-egress asymmetry — in particular the pre-ingress 
brightening — evinced in the Kepler light curve (Figure 
[3| finds no ready explanation under the dual planet hy- 
pothesis. 

Another way to induce precession would be to have 
a second planet orbit the first in a binary planet con- 
figuration (Podsialowski et al. 2010). However, we have 
verified, using the dynamical stability criteria of Mikkola 
(2008, and references therein), that for the parameters of 
KIC 12557548, two Jupiter-sized planets cannot be in a 
stable mutual orbit and avoid Roche lobe overflow. 

3.3. Occultation by Debris from a 
Disintegrating Rocky Planet 

Here we consider another scenario — that the occul- 
tations are produced by a much smaller, rocky planet, 
heated to such high temperatures that it is vaporizing. 
In evaporating away, it produces a debris field of dust 
that variably obscures up to ~1% of the light from the 
parent K star. 

A purely gaseous wind emitted by a planet will not 
have sufficient broadband opacity to affect significantly 
the transit depth in the Kepler bandpass. For example, 
the hot Jupiter HD 209458b is known to emit a gaseous 
wind (Vidal-Madjar et al. 2003; Ben-Jaffel 2008; Linsky 
et al. 2010). Although this wind increases substantially 
the optical depth in narrow atomic absorption lines, it 

5 For an assumed orbit with P = 2 X 15.685 hr, and two essen- 
tially equal eclipses per orbit, the limit on the mass of the com- 
panion would be ~11 Mj (for M p <C M s ). 
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hardly alters the transit depth in the broadband opti- 
cal. For KIC 12557548, we must consider instead that 
the planetary outflow contains particulates, e.g., dust. 
Such dust would be driven off a rocky planet so hot that 
it is vaporizing. This scenario is further developed and 
quantified in S|4] 

3.4. Other Orbital Configurations 

Could the source of the occultations be a low-mass 
eclipsing stellar binary that either is coincident ally lo- 
cated close on the sky to KIC 125 5754 8 (but see the 
limits we have placed on blends in j |2.4| ) or orbits KIC 
12557548 in a hierarchical triple system? Such a binary 
might consist of, e.g., a pair of M stars, whose strongly 
modulated light would be diluted by light from the K star 
seen in KIC 12557548. If such a binary had an orbital 
period about the K star of >10 days, it could avoid pro- 
ducing a detectable ELV in the light curve of the parent 
star. But we know of no mechanism by which an eclips- 
ing stellar binary like this, whether it is or is not in a 
hierarchical triple, can generate variable transit depths 
(e.g., orbital precession is too slow), let alone reproduce 
the transit profile shown in Figure 3 with its distinctive 
pre-ingress brightening and prolonged ingress. 

Another possibility is that one component of the low- 
mass eclipsing stellar binary is a compact object, e.g., a 
white dwarf, fed by an accretion disk (G. Marcy, personal 
communication). Light from the binary, made variable 
by occultations of the accretion disk by the donor star, 
would be diluted by light from the K star. To explain 
the variable transit depths, the accretion disk itself would 
have to vary in luminosity. Why such luminosity varia- 
tions would be restricted to those phases when the disk 
is occulted — i.e., why they would not also occur during 
the ~70% of the orbital cycle when the light curve is 
essentially flat — is not explained under this hypothesis. 

4. OBSCURATION BY DUST FROM A DISINTEGRATING 
SUPER-MERCURY 

In this section we examine in some detail how dust 
emanating from a possibly rocky planet might result in 
variable occultations of the parent star of up to 1.3% in 
depth. We discuss the size of the required dust cloud, 
composition and survival requirements for the dust par- 
ticles, and required mass loss rates and implied lifetime of 



the planet against evaporation (§[4.1 4.2 



4.3 1. A mech- 



anism of ejecting gas and dust from a planet with non- 
negligible gravity is described in §4.4| We argue that dust 
can provide a natural explanation for the time variability 
of the m ass l oss and thereby of the variable occultation 
depths (j |4.5[ ). We carry out a numerical simulation of 
the motions of dust particles after they escape from the 
planet's gravitational well that shows that a comet-like 
tail is formed, and that illustrates how t he dust might 
scatter radiation as well as absorb it (§^4.6 4.7 1. We 
go on to discuss how long such a close planetary system 
could las t under the influence of tidal drag on the or- 
bit Q4.8). Finally, we explain why occultations by dust 



have not been heretofore seen in other close-in, hot, rocky 
planets (jp~9|. 



A number of these sections are necessarily speculative 
or involve some assumptions. We offer these simply to 
start the thinking process on this subject rather than to 
provide definitive answers. 



4.1. The Underlying Planet vs. 
the Size of the Dust Cloud 

To set the stage, we assume in what follows that the 
planet underlying the outflow (hereafter "KIC 1255b") 
has an Earth-like composition, and has mass M p = 
0.1 JWq (1.8 times the mass of Mercury) and radius 
Rp — 0.5 i?©. In the absence of an outflow, our nom- 
inal "super-Mercury" would produce a transit depth of 
at most (R p /R s ) 2 « 5 x 10~ 5 , consistent with our lack of 
detections of transits in portions of the light curve. The 
planet's Hill sphere radius (distance from the planet to 
the nearest Lagrange point) is 



( M p 



\3M S 



1/3 



1 x 1Q 



0.1 M ffi 



1/3 



cm . (4) 



The Hill sphere radius is smaller than i? Q , defined as 
the characteristic size of the occulting region during the 
deepest eclipses (i.e., the size of the dust cloud). We 
adopt a nominal value for R a = 0.1 i? s ~ 5 x 10 9 cm, 
which is such that if the occulting region were circular 
and optically thick, it would generate a transit depth 
of 1%, near the maximum observed. From our experi- 
ments (not shown) simulating the occultation profile, we 
estimate that R a could be a few times larger than our 
nominal choice if the occulting region were optically thin 
and/or the transit were not equatorial. 

We have chosen a Mercury-like planet for our fiducial 
case because a low escape velocity allows us to explain 
how dust might be ejected. The physics of how grains 
can be launched successfully from a super-Mercury is ex- 
plained in |4.4| In contrast, the difficulties of eje cting 
dust from hot rocky super-Earths are discussed in £ 4.9 



4.2. Dust Composition and Survival 
Dayside planet temperatures peak at approximately 



! cff,p 




2100 K, 



(5) 



neglecting, among other effects, transport of heat from 
the substellar point by the atmosphere or by a magma 
ocean (see Leger et al. 2011 for an estimate of the location 
of the shore of the magma ocean for the hot super-Earth 
CoRoT-7b). Such nominal temperatures are hot enough 
to vaporize silicates and a variety of solid oxides, and to 
form a high-Z atmosphere. 

This particular model requires that the atmosphere 
contain refractory grains which can survive for long 
enough to produce an observable occultation. We can 
think of at least two ways of loading the atmosphere with 
dust. The first is by condensation of dust grains from 
a metal-rich vapor into clouds akin to those in brown 
dwarf or giant planet atmospheres; see, e.g., the model- 
ing of forsterite clouds in Cooper et al. (2003) and Mad- 
husudhan et al. (2011). Condensation can occur where 
gas cools, at altitude or near the day-night terminator 
(Schaefer & Fegley 2009; Schaefer et al. 2011; Castan & 
Menou 2011). Whether such clouds can contain eno ugh 
condensed mass remains to be seen — we will see in |4.3| 
and j j4.4| that the planetary wind must contain at least 
as much mass in gas as in solids (or liquid droplets). 

A second way to inject dust into the atmosphere is 
through explosive volcanism. Tidal locking of hot rocky 
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planets may produce a pattern of mantle convection in 
which volcanism is foc used near the substellar point, 
where a thermal wind (^4.4) is best able to entrain vol- 
canic ash (Gelman, Elkins- Tanton, & Seager 2011; van 
Summeren, Conrad, & Gaidos 2011). Weakening and 
thinning of the planetary lithosphere by sublimation may 
render the surface particularly vulnerable to volcanic ac- 
tivity. Eruptions beyond the shore of the magma ocean 
may be supplemented by violent bursting of bubbles at 
the magma ocean surface: on Earth, such eruptions 
can eject magma spray up to ^400 m/s (Taddeucci et 
al. 2012). Io's plumes are a Solar System analog for vol- 
canic injection of dust into a tenuous atmosphere; dust 
grain radii range from 0.03-0.12 /im (Jessup & Spencer 
2012). 

For completeness, we note that there might also be a 
third way to pepper the atmosphere with dust: as the 
still solid portions of the planetary surface sublime, the 
resulting vapor may entrain grains of more refractory 
composition. 

The upper mantle of the Earth is composed over- 
whelmingly of pyroxene ([MgjFeJSiOs) and olivine 
([Mg,Fe] 2 Si04). We might expect whatever grains to be 
present in the atmosphere of an evaporating sub-Earth 
to also consist of these minerals, at least in partj^] The 
observations of KIC 12557548 demand that grains not 
sublimate before they can cover up to ~1% of the face 
of the star. This requirement can be met by micron- 
sized grains of pyroxene, but not of olivine. In a study 
of grain survival times in the tails of sun-grazing comets, 
Kimura et al. (2002; see their Figure 4) calculated that 
the sublimation lifetime of an amorphous pyroxene grain 
having a radius s ~ 0.2 /im is i su b ~ 3 x 10 4 s at a dis- 
tance of 7.5 Rq from the Sun (equivalent to a distance 
of 2.8 Rq from the K star in KIC 12557548, neglecting 
the difference in the shape of the stellar spectrum heat- 
ing the grain). Crystalline pyroxene survives for still 
longer, a consequence of its lower absorptivity at optical 
wavelengths and hence cooler temperature. By contrast, 
micron-sized grains of olivine ([Mg,Fe]2Si04) vaporize in 
mere seconds FH 

A sublimation lifetime of t S ub > 3x 10 4 s, as is afforded 
by micron-sized pyroxene particles, is long enough for 
grains to travel the length of the occulting region. This 
travel time is on the order of 



sphere: 



''travel 



Ro 

V 



2 x 10 4 s 



(6) 



where the first term accounts for the time to travel from 
the planetary surface to the Hill sphere, and the second 
term accounts for the time to travel out to R . For the 
first leg of the journey, we have adopted a flow speed 
equal to the sound speed c s in the planet's high-Z atmo- 



16 Another possibility involves grains of alumina (a.k.a. corun- 
dum; AI2O3), which seems likely if the planetary surface has 
been distilled to the point where it is composed predominantly 
of calcium-aluminum oxides (see Figure 5 of Leger et al. 2011). 

17 We have been unable to identify a simple reason why the 
sublimation times of olivine grains are orders of magnitude shorter 
than those of pyroxene grains; see references in Kimura et al. (2002) 
for the original laboratory data. 



/ kT 



0.7 



/ T 



V2000K 



1/2 / 30 \ 1/2 



km/s (7) 



where k is Boltzmann's constant, p is the mean molecular 
weight, and mjj is the mass of hydrogen. The sound 
speed is appropriate to use if the grains are carried off 
the planetary sur face by a thermal Parker-type wind, as 
discussed in §4.4| For the second leg of the journey, we 
adopt v ~ 10 km/s. This is a characteristic grain velocity 
once the wind is rarefied enough that dust decouples from 
gas, and is given by the dyna mics dictated by gravity and 
radiation pressure (see ^4.6). 



4.3. Mass Loss Rate and Evaporation Timescale 
The mass loss rate in dust grains may be estimated as 



M d ~ Qp d vRl 



(8) 



where pd is the mass density in dust grains at a distance 
R from the planet is the outflow velocity (on scales 
of R a ), and f2 is the solid angle subtended by the flow 
as measured from the planet center. We can relate pd to 
the optical depth through the dust flow, integrated along 
the line-of-sight path of length 1Z ~ R a : 



T ~ PdR-Kd 

where the dust opacity Kd ~ 



PdRo 
PbS 



(9) 



;2 /(PbS 3 ) ~ l/(pbs) for 
grains of radius s and bulk density pb ~ 3g cm -3 . We 
can also relate r to the eclipse depth /, assuming r < 1 
(we will discuss why t he flow may be optically thin, or 
marginally so, in S4.5): 



/ 



t_ri 



0.01 



(t) 



R0/R1 
0.1 



(10) 



Combining the above three relations, we have 



A v gft#» a 2 x 1Qll fd 



Ro 

f 
0.01 



lM e /Gyr 



Pb 



1/ \3g/crn / \0.2pm 
10 km/s ^ \R /R S ) \0.65R Q J 8/& 



(11) 

where fi is normalized to 1 to account for the likelihood 
that the area from which dust is launched ma y pr imarily 
be in the substellar region. If, as noted in |4.1| t < 1 
and/or the transit is not equatorial, then R a may be up 
to a few times larger than 0.1i? s , and the above estimate 
of Md would need to be reduced somewhat. 

Gas (metal-rich vapor) adds to the total mass loss rate 
by a factor of £. If we assume that the mass loss mech- 
anism principally produces a gaseous outflow, and that 
drag accelerates the dust grains, then it seems reasonable 
to assume that the gas density must be at least as large 
as the dust density: £ > 2. Otherwise the backreaction 
drag on gas by dust would slow the gas down and pre- 
vent it from escaping. We will revisit the factor of £ more 

18 For this back-of-the-envelope estimate, we neglect the differ- 
ence between the planet center and the center of the dust cloud. 
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quantitatively in §4.4| for now, we note that £ ~ 2 falls 
squarely in the range of values inferred for cometary out- 
flows of dust mixed with gas (Delsemme 1982; Fernandez 
2005). Then the evaporation timescale for our nominal 
M p = 0.1 M e planet is 



t. 



M„ 



cvap 



0.2 



0.3 



Gyr (12) 



where T is the fraction of time the planet spends losing 
mass at the maximum rate (set by / = 0.01). Our esti- 
mate of tevap is nominally shorter than the likely age of 
the star, but not by an implausibly large factor. 

4.4. Ejecting Dust Grains Via a Thermal Wind 

In this section, we explore the possibility that dust 
is ejected from KIC 1255b in a manner that resembles 
dust ejection from comets (for a pedagogical review of 
comets, see Rauer 2007). Gas sublimates and then flows 
off a comet nucleus at essentially the sound speed, which 
can be dozens of times greater than the escape veloc- 
ity from the nucleus. As gas atoms stream toward the 
cometopause — the surface of pressure balance between 
cometary gas and the Solar wind — they are photoion- 
ized by stellar ultraviolet photons. At the cometopause, 
cometary ions are rapidly picked up by the magnetic field 
of the Solar wind, and accelerated up to the Solar wind 
speed of hundreds of km/s. This explains the well-known 
observation that the ionized gas tails of comets are swept 
anti-Sunward, in the direction of the Solar wind. Suffi- 
ciently small dust particles, on the other hand, escape 
the nucleus by virtue of the hydrodynamic drag exerted 
by the escaping gas. 

An important difference between a comet and KIC 
1255b is the much higher surface gravity of our nomi- 
nal super-Mercury. The surface escape velocity 



M„ 



0.1 M a 



1/2 



0.5i?g 



1/2 



km s 



(13) 



is ~7 times higher than the local sound speed c s . At first 
glance, then, it seems unlikely that gas can be driven off 
the surface of the planet. But in Parker-type winds (see, 
e.g., Lamers & Cassinelli 1999) the base of the wind can 
be characterized by just such values of f e sc,suri /c s ■ Start- 
ing at their base, Parker-type winds accelerate by gas 
pressure over long distances — several planetary radii — 
before the bulk wind speed v formally exceeds the lo- 
cal escape velocity. For example, u e sc,surf/c s ~ 5 in the 
Solar corona, which is the seat of the Solar wind (e.g., 
Lemaire 2011). Another example is provided by the ther- 
mal winds from hot Jupiters, where i> e sc,surf /c« ~ 5 at the 
base of the planetary wind where stellar UV photons are 
absorbed and heat the upper atmosphere (e.g., Murray- 
Clay, Chiang, & Murray 2009). Note that both these 
examples represent flows that are continuously heated 
and energized along their extents — by photoionization in 
the case of the hot Jupiter wind, and by conduction and 
wave heating in the case of the Solar wind (e.g., Marsch, 
Axford, & McKenzie 2003). Thus these winds can remain 
nearly in the hydrodynamic regime and enjoy large mass 
fluxes, even at the large Jeans parameters (Ao = v1 ac /2c 2 s ) 
characterizing their base radii -Rbasc 

19 To quantify this statement, note that a hot Jupiter wind 



The thermal winds noted above suggest by analogy 
that the atmosphere of KIC 1255b can be gradually accel- 
erated by gas pressure to speeds that eventually exceed 
the local escape velocity. Gas can be continuously heated 
along the extent of the wind by collisions with grains, 
which in turn are heated by starlight; thus the wind 
should maintain a nearly constant temperature > 2000 
K. At distances of Rsm ~ 3i? p , near the first and second 
Lagrange points where the effective gravitational accel- 
eration is nearly zero, we anticipate that the outflow will 
have attained the sonic point, so that the wind velocity 
will be of order c s on these scales (see, e.g., Figure 9 of 
Murray-Clay et al. 2009). This explains the first term in 
equation [6j Of course, this sketch needs to be confirmed 
with a first-principles calculation of the mass loss rate as 
a function of planet mass. 

How dense and fast must the gaseous outflow be to lift 
grains off the surface of the planet? The aerodynamic 
drag force on a grain must exceed the force of planetary 
gravity: 

(14) 



PgV, 



<e s s 2 <; PbS 3 g 



where g ~ 400 cm s~ 2 is the surface gravitational ac- 
celeration, and p g and v g are the gas density and flow 
velocity, respectively. For the drag force on the left-hand 
side we have used the Epstein (free molecular) drag law, 
appropriate for grains whose sizes are less than the colli- 
sional mean free path between gas molecules, and for flow 
velocities v g < c s (Epstein 1924). Both of these condi- 
tions can be verified a posteriori to hold throughout the 
flow. The atmospheric density just above the ocean melt 
is p g r~j P/Cg, where the surface pressure P ~ 10 ph&r 
at T = 2100 K, if the ocean is composed predominantly 
of Mg, Si, and O, with zero contribution from Na (see 
Figures 5 and 6 from Leger et al. 2011, and consider va- 
porized fractions 0.85 > F vap > 0.02; see also Schaefer 
& Fegley 2009). Such an atmosphere is composed pri- 
marily of Mg, SiO, O, and O2 gas, according to Figure 
3 of Schaefer & Fegley (2009). Then the vertical wind 
velocity required to keep grains aloft against gravity is 



10 /xbar 
P 



0.2 pm 



m/s 



(15) 



which seems easy to satisfy. If the surface has been va- 
porized (distilled) to the point where it is composed of 
Al, Ca, and O (F vap > 0.85), then P ~ 0.1 pb&r and 
v g > 200 m/s, which still seems possible (at least hori- 
zontal wind velocities are thought to have near-sonic ve- 
locities; Castan & Menou 2011). 

The lift-off condition ( 14 1 implies a mass loss rate in 
gas of 



Mg ~ npgVgR 2 p > 4 x 10 1 



0.2 (im 



g/s. (16) 



(as computed, e.g., in the standard model of Murray-Clay et al. 
2009) is characterized by < t'/4 , o 1 o ~ 0.05 at Ao ~ 13.5, in the no- 
tation of Volkov et al. (2011). Here "t'/'t'o.O is the dimensionless 
mass loss rate, and the "0" base radius is taken to be the surface 
which presents unit optical depth to UV photons. For compari- 
son, Volkov et al. (2011), who neglect energy deposition at radii 
> ^basei predict &/$ofi ~ 4 X 10~ 5 at Ao = 13.5 (see their Figure 
3). Thus the continuous energy deposition enjoyed by a hot Jupiter 
wind increases the mass loss rate above what it would be without 
continuous energy deposition by about three orders of magnitude. 
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The fact that the minimum value of M g is comparable to 
Md as estimated in (11) supports our claim in £4.3 that 
the gas-to-dust ratio may be of order unity, i.e., £ ~ 2. 

It is remarkable that at a maximum eclipse depth / = 
0.01, our order-of-magnitude estimate for the total mass 
loss rate £Md ~ 4 x 10 11 g/s from our super-Mercury 
either matches or exceeds, by up to an order of magni- 
tude, estimates of mass loss rates — mostly in hydrogen 
gas — from hot Jupiters like HD 209458b (Vidal-Madjar 
et al. 2003; Yelle 2004; Garcia-Munoz 2007; Murray- 
Clay et al. 2009; Linsky et al. 2010) and HD 189733b 
(Lecavelier des Etangs et al. 2010). Hot Jupiters are 
thought to lose mass by photoevaporation, whereby stel- 
lar X-ray and ultraviolet (XUV) radiation photoionizcs, 
heats, and drives a thermal wind off the uppermost layers 
of a planetary atmosphere. Present-day mass loss rates 
of hot Jupiters are limited by the XUV luminosities of 
main-sequence solar-type stars, of order L X uv,s ~ 10~ 5 - 
1O _6 L0. In the case of KIC 1255b, mass loss may also be 
driven by a thermal wind — but one which is powered by 
the bolometric stellar luminosity, L s ~ O.IALq, which 
heats the entire planetary atmosphere to temperatures 
T ~ 2000 K. The power required to drive our inferred 
mass loss rate 



L 



required ~ 



GM p £M d 



>6 x 10 22 ( | 



2 / \0.2fim 



(j) er s/ s 



(17) 

(G is the gravitational constant) is easily supplied by the 
total power intercepted in stellar radiation, 

^intercepted = '4~ 4 X 10 26 erg/s . (18) 

4.5. Time Variability of Mass Loss and 
Limiting Mass Loss Rate 

The variability observed in the obscuration depths im- 
plies substantial changes in Md (or in other physical 
properties of the outflowing dust clouds) over timescales 
< P or b = 15.7 hr. Although we cannot be quantitative, 
such rapid variability seems possible to accommodate 
within the framework of an evaporating planet. Again, a 
comparison with comets is instructive. As in a cometary 
nucleus, inhomogeneities or localized weaknesses in the 
planet's surface ("vents"; Jewitt 1996; Hsieh et al. 2010) 
may cause the planet to erupt occasionally in bursts of 
dust. Whatever gas is ejected from the planet will be 
photoionized by the stellar radiation field and, like the 
gas tails of comets, interact with the magnetized stellar 
wind. To the extent that dust can be entrained in plasma 
along the planet's magnetotail, plasma instabilities may 
introduce further time variability. 

The mass loss rate Md may be limited to ^2 x 10 11 g 



the value in equation (11) appropriate to an opti 



cal depth r ~ 1 (assuming aTTother parameters in that 
equation are fixed at their nominal values). This upper 
limit exists because if the flow became too optically thick, 
light from the star would not reach the planet's surface to 
heat it. Thus the mass loss rate may be self- limiting in a 
way that keeps the outflow from becoming too optically 



thick. One can imagine a limit cycle which alternates 
between a low-obscuration phase during which starlight 
strikes the planet's surface directly and dust proceeds 
to load the planet's atmosphere, and a high-obscuration 
phase during which dust removal dominates production. 

Wind mass loss rates should depend on outgassing 
rates from the planetary surface. Since saturation va- 
por pressures depend exponentially on temperature, it 
seems easy to imagine that modest changes in surface 
temperature due to changes in the optical depth of the 
wind will drive large changes in Md- As noted above, a 
first-principles calculation of M d along the lines of the 
externally irradiated solution for hot Jupiter winds (e.g., 
Murray-Clay et al. 2009) would be welcome. 

4.6. Dust Stream Structure 

Once dust escapes from the planet, and the gas density 
becomes too low to exert significant drag, the dust will 
be swept back into a comet-like tail by the combination 
of radiation pressure and Coriolis forces. The ratio of 
the stellar radiation pressure force on a grain to that of 
stellar gravity equals 



= 



*rad 



3Q pr -L 



16irGcM s p h s 
L s \ / O.7M 
0.14 L J \ M s 



:0.15 



0.2 peax 
s 



(19) 



where Q pl - is the radiation pressure efficiency factor (see, 
e.g., Burns et al. 1979) and c is the speed of light. For 
the numerical evaluation in ( 19 ), we have used the results 



of Kimura et al. (2002, their Figure 2) for a spherical 
pyroxene particle, scaled for the lu mino sity and mass of 
our K star. The numerical result in ( 19 ) is valid only for 



> 0.2 ^m (as s decreases below 0.2 yim., Q pr decreases 
and radiation pressure is increasingly negligible; in any 
case such grains are unimportant because they sublimate 
too quickly). 

We performed a simple calculation of what a dust tail 
might look like. Since the precise details of how dust 
is launched from the surface are presently unclear — they 
are subject, e.g., to horizontal winds on the planet that 
can distribute material from the substellar point to the 
nightside — the calculation simply assumes a spherically 
symmetric outflow of dust from the planet at the sur- 
face escape speed. Our fiducial planet parameters of 
Rp = 0.5 R m and M p = 0.1M 9 , /3 = 0.15, and a sub- 
limation time for the dust of i S ub = 10 4 s were assumed. 
The computations were performed in a frame of reference 
rotating with the planet's orbit, and followed 30,000 dust 
particles from launch to sublimation. 

The results of this simulation are shown in Fig.[6| A 
long dust tail is apparent because the color scaling covers 
a large dynamic range of column densities. The effective 
length of the tail is not actually that long, e.g., the col- 
umn density of dust particles is quite small at distances 
comparable to R s . This was quantified by integrating 
the total number of dust grains within ± 1/8 of the stel- 
lar radius in height (z) as a function of distance along 
the dust tail (y). The result is shown in Fig. [7] Given 
the presently available information, only relative values 
of the optical depth can be shown, but the curve serves 
to show that the projected density of dust grains falls 
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Figure 6. Simulations of dust flows from a 0.1 M® planet of 0.5 
driven by radiation pressure from the parent star. The spatial 
scale is in units of the semimajor axis of the planet. The color 
scale is proportional to the square root of the projected density in 
order to enhance the dynamic range. The red circle represents the 
parent star. The top panel shows the dust density as viewed from 
the orbital pole; the bottom panel shows the dust density as seen 
by an equatorial observer at mid transit. 



off sharply with distance from the planet. Thus, the op- 
tical depth could be somewhat greater than unity near 
the projected planet disk, and the net obscuration of the 
parent star could take on values up to the maximum level 
of ~1.3% that is observed. 

In summary, the dust stream geometry seems able to 
accommodate the occultation durations and depths. It 
might even be able to explain the curious ingress-egress 
asymmetry in the occultation profile, as we argue in the 
next section. 
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Figure 7. Relative optical depth of the simulated dust tail as a 
function of distance across the parent star in units of the semimajor 
axis. This result is derived from the simulated data used to produce 
Fig. [6] The vertical scale is arbitrary since otherwise we would need 
to know the absolute value of the optical depth near the surface of 
the planet. 



4.7. Scattering 

The putative dust cloud likely attenuates more by scat- 
tering photons out of our line of sight than by actual 
absorption (this is due to the small imaginary part of 
the index of refraction in the visible; see, e.g., Denevi et 
al. 2007). This opens the possibility that scattering of 
starlight into the line of sight would produce, at certain 
orbital phases, an excess of light above the nominal out- 
of-occultation region levels. An excess would be more 
likely given the presence of larger (s > 1 /im) grains that 
scatter primarily in the forward direction. 

We can estimate the ratio of fluxes scattered into, ver- 
sus out of, our line of sight. Consider an idealized, small, 
spherical scattering cloud of radius Rd that scatters a 
fraction A (for "albedo") of all the stellar radiation in- 
cident upon it (i.e., the cloud presents a total scattering 
cross section of A ■ irlQ)- The incident radiation is scat- 
tered uniformly into a solid angle £l sca t centered about 
the direction of the incident beam. In the case that the 
cloud sits directly in front of the target star, the ratio of 
fluxes scattered into, versus out of, the line of sight is 



Xin 
Xout 



i6 a 



(20) 



independent of A and Rd- The ratio ranges from ~1.6% 
for isotropic scattering, to as high as ~20% for forward 
scattering into 1 steradian. If the dust cloud is blocking 
the stellar disk, then attenuation is dominant (except in 
the case of extreme forward scattering) . However, if the 
dust cloud is not blocking the stellar disk, then there may 
be a small excess contribution to the system light. This 
may explain the small bump in flux just before ingress, 
and the small depression in the flux just after egress (see 
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Fig. [3]). In this interpretation, the ingress-egress asym- 
metry results from the asymmetric distribution of dust 
in the planetary outflow. Just prior to ingress, light is 
forward-scattered off the dusty head of the "comet" into 
our line of sight, producing the small bump in flux. Light 
may continue to be forward scattered into the line of sight 
during the main occultation and even after egress. How- 
ever, at those phases when part of the tail occults the 
star, the scattering of starlight out of the line of sight is 
likely to dominate and produce a net depression in the 
observed intensity. 

4.8. Tidal Inspiral Lifetime 

We now consider how long such a planet as the one we 
are proposing might last in the face of tidal decay due to 
interactions with the parent star. The timescale of tidal 
decay is approximately 



T tidal — Q 



2M S 
9M„ 



P, 



orb 



2tt 



10 



10 



Q_ 
10 6 



yr 



(21) 



(see, e.g., Gu, Lin, & Bodenheimer 2003) where Q is the 
usual dimcnsionless measure of dissipation in the parent 
star, and we have assumed that the parent star is rotating 
with an angular frequency much lower than the planet's 
orbital frequency, and that the orbital eccentricity is zero. 
Thus, even for a pessimistically small value of Q, the 
lifetime against tidal decay is more than adequately long. 

4.9. Differences Between KIC 12557548 and Hot 
Super-Earths 

If our proposed scenario for KIC 1255b is correct, 
why don't other hot rocky planets also evince occult- 
ing dust clouds? The class of close-in rocky planets in- 
cludes CoRoT-7b (Leger et al. 2009; Leger et al. 2011), 
55 Cnc e (Winn et al. 2011), and Kepler-lOb (Batalha 
et al. 2011b)Pj All are argued to have significant rock 
components, although an admixture of water, hydrogen, 
and/or other light elements seems required to explain 55 
Cnc e. 

The primary parameter distinguishing the super- 
Mercury that we are imagining orbits KIC 12557548 and 
the other hot rocky planets is clearly the planet mass. 
All the other planets are "super-Earths" having masses 
M p ~ 10 M® and radii R p ~ 2i? (mass-radius rela- 
tions for super-Earths, may be found in, e.g., Valencia et 
al. 2007 or Figure 3 of Winn et al. 2011). Super-Earths 
have surface escape velocities that, at w CSCiP ~ 25 km/s, 
are almost certainly too large for gas that is heated to 
temperatures T ~ 2000-3000 K to escape at substantial 
rates via Parker- type winds. The efficiency of driving a 
thermal wind drops rapidly with increasing v csc ^ p . The 
local hydrostatic scale height of gas (valid for altitudes 
h < RJ 



2c 



^scl 



-Rr, 



2 x 10" 



/ 10Ma 



Rp 
2R a 



R 



p 
22) 



20 Kepler-9d (Batalha et al. 2011a) may also belong to this class 
based on its size and orbital distance, but, unfortunately, its mass 
is not yet known, and therefore we cannot be sure that it is a purely 
rocky planet. 



decreases as l/v^ sc . For super-Earths, the base of a 
Parker-type wind (the base is roughly identified as that 
location where the local value of v esc /c s ~ 6) would lie 
several R p away from the planetary surface, while for a 
super-Mercury the base would lie close to the surface. 
Because the gas density must decrease more or less ac- 
cording to p g cx exp (— h/h sc \) from the surface to the 
base, the gas density at the base on a super-Earth would 
be many orders of magnitude lower than on a super- 
Mercury. 

Alternate mechanisms of mass loss from super-Earths, 
other than a Parker-type wind, include: (i) XUV pho- 
toevaporation (by analogy with hot Jupiters) , (ii) stellar 
wind drag, and (iii) stellar radiation pressure. All fail 
to explain the inferred mass loss rates for KIC 1255b. 
The intercepted power from the stellar XUV luminosity 
is at least a fact or of 10 too low compared to £ req uired 
(see equation 17 1, even for optimistically high values of 
ixnv.B = lO^^L^ and of the efficiency e — 0.1 with 



J XUV,s 

which XUV 
cal work 



nhoton energy is converted to mechani- 
Stellar wind drag suffers the same prob- 
lem that XUV photoevaporation does: the power in- 
tercepted by the planet from the stellar wind is insuf- 
ficient to supply ^required- An easy way to see this is to 
note that the mechanical luminosity of the Solar wind is 
iwi„d, s ~ (l/2)M t£. ndiS ~ (1/2) • (2 x lO" 14 M /yr) • 

23J24 



The inter- 
to account 



(400 km/s) 2 - 3 x 1O- 7 L < L X uv, s _ 
cepted power is still too low even were we 
for a planetary magnetosphere, which would increase the 
cross section presented by the planet by a factor < 10r^ 
Finally, based on the calculations of Kimura et al. (2002J7 
the force of stellar radiation pressure acting on a grain 
is at most ^30% of the force of planetary gravity on 
a super-E arth, for conditions in KIC 12557548 (but see 
§S 4.6-4.71. Moreover, even if the force of radiation pres 



sure on a grain were somewhat larger than that of plane- 
tary gravity, the grain would be dragged by ambient gas, 
and might never accelerate to escape velocity. 

For all these reasons, super-Earths cannot launch dusty 
outflows like those we infer for KIC 12557548 — whereas 
sub-Earths can. There is still another reason why dust 
clouds are especially observable for KIC 12557548 and 
not in the other observed systems. Even if dust grains 
could be driven off the hot super-Earths detected to date, 
they would be too hot; grains in these other systems 
would sublimate too quickly to travel far from their par- 
ent planets, and would fail to block as much light from 
their host stars. If we use equation (J5]) to estimate the 
temperature of a grain at the location of each super- 

21 XUV luminosities of K dwarfs are typically 3 X greater than 
those of G dwarfs (Hodgkin & Pyle 1994; Lecavelier des Etangs 
2007). 

22 The efficiency e is lowered by ionization and radiative losses; 
radiative losses may be especially severe in the case of KIC 1255b 
because of cooling by dust, a problem that does not afflict hot 
Jupiter winds. 

23 Mass loss rates of K stars are observed to be comparable to 
that of the Sun (Wood et al. 2005; Cranmer 2008). 

24 At the orbital distance of KIC 1255b, the energy carried by 
the stellar wind is mostly magnetic and thermal, not kinetic. But 
after accounting for all three forms of energy using empirical mea- 
surements of the Solar wind at a = 2.8Rq, we arrived at roughly 
the same total power output, L w i n d. s ~ 10 — 6 Lq. 

25 If KIC 1255b has a surface dipole field of 1 G, then its mag- 
netospheric radius is ~2i?. p . 
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Earth, we find: 

• Kepler- 10b : 3000 K 

• 55 Cnc e : 2800 K 

• CoRoT-7b : 2500 K 

• KIC 1255b : 2100 K 

Thus KIC 1255b — orbiting the coolest star — is distin- 
guished in having the coolest grains. Sublimation life- 
times t su b are exponentially sensitive to temperature. 
For example, whereas we have estimated t su b ~ 3 x 10 4 s 
for 0.2-/zm-sized pyroxene grains emitted by KIC 1255b, 
these same grains emitted by CoRoT-7b would have 
t S ub ~ 10 2 s (Kimura et al. 2002) — considerably shorter 
than the travel time t * rm ^ \ ~ 2x 10 4 s required to produce 
a ~1% occultation (j |4.2| . 

Thus KIC 1255b apparently occupies a sweet spot for 
the production of occulting dust clouds: the system is 
hot enough for the surface to vaporize and for a wind to 
be launched)^] but cool enough that dust particles can 
be formed from certain minerals not quite hot enough to 
vaporize (or to form by condensation upon cooling of the 
vapor after it flows to a different location), and to allow 
these dust particles to travel a significant fraction of the 
host stellar radius. 

5. SUMMARY AND OUTLOOK 

We have reported on occultations of KIC 12557548 
that recur with a period of 15.6854 hours, have depths 
that range from 1.3% to < 0.15%, and that vary on 
timescales comparable to, if not also shorter than, the 
occultation period. A search for periodic modulation in 
the depths of the occultations did not reveal anything 
definitive. From various lines of evidence, most notably 
a newly acquired optical spectrum, the target star ap- 
pears to be a mid-K dwarf. 

We briefly considered a scenario wherein the changing 
transit depths are due to a dual giant planet system (ei- 
ther on separate orbits or in a binary pair) in which one 
of the planets undergoes grazing transits of the target K 
star. The orbital plane of the transiting planet would 
have to precess in such a way as to effect the variable oc- 
cultation depths. We find, however, that a binary planet 
configuration is very likely to be dynamically unstable, 
even if the two planets are themselves essentially in a 
contact configuration. In the case of two planets on sep- 
arate orbits, it is not clear how orbital precession could 
explain either the remarkable changes from orbit to or- 
bit or the long intervals of ~10 and 20 days when the 
occultations mostly disappear. We also briefly consid- 
ered an eclipsing binary that is orbiting KIC 12557548 
in a hierarchical triple configuration, but were unable to 
explain the basic properties of the observed light curve 
in the context of such a scenario. 

We come down in favor of a scenario in which dust is 
aerodynamically dragged off an orbiting rocky planet — 
possibly one not much larger than Mercury in size — by 
a thermal Parker-type wind composed of metal atoms 
sublimated off the planet's surface at a temperature of 

26 The three sub-Earths orbiting KOI-961 (Muirhead et 
al. 2012) are probably too cold to satisfy this condition. 



^2000 K. The atmosphere may be loaded with dust ei- 
ther through condensation into clouds, or through explo- 
sive volcanic activity. The mass loss rate must sometimes 
be as high as a few times 10 11 g/s, or about 1 M e /Gyr, 
and may comprise by mass roughly equal parts gas and 
dust. Our fiducial mass and radius for the underlying 
planet are 0.1 Mq and 0.5 i?©, respectively. The corre- 
sponding lifetime against evaporation is ~0.2 Gyr. This 
is probably less than the age of the star, but is not alarm- 
ingly short. Because the mass loss rates are inferred from 
the observations and as such are largely fixed, planets 
much smaller in mass than our fiducial super-Mercury 
have implausibly short evaporation lifetimes, while plan- 
ets much larger in mass and thus in surface gravity do 
not seem capable of sustaining the required mass loss 
rates. 

After leaving the planet's gravitational well, the dusty 
gas will soon become rarified to the point that the dust 
and gas trajectories diverge. We have simulated how 
dust grains, decoupled from gas, would be shaped into 
a comet-like tail by the Coriolis force and stellar radia- 
tion pressure. The head-tail asymmetry of the dust cloud 
promises to explain the observed ingress-egress asymme- 
try in the occultation profile. Pyroxene grains, each a 
fraction of a micron in size, are a good candidate for 
the obscuring grains because of their relatively long life- 
times against sublimation — alumina is also a promising 
high-temperature mineral to consider. Time variability 
of the outflow's dust density as reflected in time vari- 
ability of the occultation depths could be caused by any 
number of factors: irregularities and uneven outgassing 
rates on the planet's possibly volcanic surface; a feed- 
back loop between mass loss and the stellar radiation 
flux received at the planet's surface that causes both to 
fluctuate with time; and instabilities resulting from the 
interaction of the ionized planetary wind with the magne- 
tized stellar wind. The analogy between our evaporating 
super-Mercury and comets may not be a bad one, but 
much theoretical work remains to be done to put this 
scenario on a more quantitative and secure footing. 

5.1. Observational Prospects 

Obviously, short-cadence observations with Kepler 
would be helpful in clarifying the issue of whether there 
are shape changes in individual occultations vs. only 
changes in the overall depths. Some shape changes would 
be expected, as the timescale for the wind to "refresh" is 
given by the grain travel time travel ~ 5 hr (equation Rjl), 
which is not much longer than the eclipse duration of 1.5 
hr. The target star, at K p = 15.7 magnitudes, is not very 
bright, but each 1-minute sample with K epler would yield 
1 a ~ 0.15% counting statistics. Given that the depths of 
the deepest occultations are 1.3%, this would provide a 
signal-to-noise of 8:1 for individual samples. The ingress 
and egress should be resolved in short-cadence data and 
the shapes should give important clues as to the nature 
of the system. 

Observations with a large optical/infrared telescope, 
either from the ground or from space, would also prove 
helpful. Measurement of the wavelength dependence of 
the occultation depth could confirm our hypothesis that 
sub-micron-sized dust is responsible. The detection of 
solid-state absorption features (e.g., the 10 /xm silicate 
band) would also clinch the case for dust. Deep imaging 
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would have the added benefit ol further ruling out blends 
with faint, variable interlopers. 

Searching for extra in-transit absorption by metals — 
e.g., photoionized Mg, O, Si, Ca, and Fe — through ob- 
servations of spectral lines could lead to a direct demon- 
stration that the underlying planet is comprised largely 
of heavy elements, and is evaporating (cf. Linsky et 
al. 2010). Our rough calculations indicate that in some 
lines (e.g., the Mg II 2800 A doublet), the gaseous out- 
flow may be so optically thick that its occulting size is 
limited only by the confining pressure of the incident 
stellar wind. We estimate the corresponding maximum 
flux decrement in any spectral lines to be on the order of 
-5%. 

If there indeed turns out to be a super-Mercury orbit- 
ing KIC 12557548, where the occultations are due to dust 
carried by outflowing sublimated gas, the planet would 
be among the smallest-mass bodies ever to have been (in- 
directly) detected. It would be the first extrasolar planet 
shown to be geologically active and to be disintegrating 
via the loss of high-Z material. 
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